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Discrepancy between periodic orbit theory and numerical calculation ol a modified Kepler problem 
is cleared by a quantum mechanical calculation. The diagonal approximation already gives a good fit 
for the numerical calculation. A better result yet is gained by considering the coherent interference 
between the classical periodic orbits and the Balian-Bloch term. This approach produces improved 
results for the rectangular billiards as well. 
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I. INTRODUCTION 

In a companion paper, [l[ we investigated level correla- 
tions in a modified Kepler (MK) problem and discovered 
two important effects: quantum-Hall-like jumps in the 
saturation level rigidity as a function of running energy 
and persistent oscillations of the level number variance 
as a function of the interval width. Unlike rectangular 
billiards (RB), where superposition of a number of 
incommensurate harmonics led to a finite lower limit of 
the variance, in the modified Kepler problem the level 
variance can achieve near zero values. 

In other words, variance can decrease with the increase 
of the interval width due to the oscillatory behavior and, 
moreover, at certain values of the interval width the num- 
ber of levels is nearly identical for each of the members of 
the statistical ensemble. This is a highly unusual circum- 
stance from a statistical viewpoint and is due to specifics 
of the strong scmiclassical spectral correlations in classi- 
cally integrable systems. 

Previously, we used a periodic orbit theory [l[ to in- 
vestigate level correlations. However, an important class 
of periodic obits - radial oscillations for a fixed angle - 
were omitted. Here we utilized a quantum mechanical 
approach that confirms a corrected periodic orbits ap- 
proach. Both are in complete agreement now and are in 
excellent agreement with the numerical simulation. 

The key ingredient of our numerical work is the use 
of parametric averaging (for instance, the aspect ratio 
of equal area rectangles for RB; for MK see discussion 
in appendix). This allows us to do statistical sampling 
for a fixed running energy in the semiclassical spectrum 
and should be contrasted to the traditional technique of 
spectral averaging in which the key effects discussed in 
this work are washed out. 



II. LEVEL CORRELATION FUNCTION AND 
LEVEL NUMBER VARIANCE 

A. average spectral staircase 

The eigen-energy of Modified Kepler problem E = 
2puj + I 2 , where oj = ^/2f3 [J. The average cumulative 



number of levels is 
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which is the scaling relation we use in the numerical cal- 
culation. But in the theoretical calculation, only the 
first term is considered. We define the scaled energy 
e = (,yV{E)) « E 3 / 2 /3uj. 

The average spectral staircase of rectangular billiard 
(RB) is given in Eq. (5) of Ref. 0. 



B. variance 

1. level fluctuation 

The level fluctuation Q 

Sp = p-(p)=Sp^+Sp^, (2) 

where Sp^ is the Balian-Bloch term explained in ap- 
pendix and Sp^ is given by 
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where the second term is due to radial periodic orbits 
and the third term due to circular periodic orbits. The 
level density fluctuation equation is the starting point to 
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calculate spectral statistics. For RB, 

5 P W = (l/^TT 3 ) 1 / 4 
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where Mi = Mia 1 / 4 and M 2 = M 2 a" 1 / 4 . 



2. Zet>eZ correlation function 

The correlation function of levels is defined as 
K(e;E) = (Sp(e - E/2)5p(e + E/2)), where the ensem- 
ble average (•) is defined as parametric averaging over an 
ensemble of w. 

Using diagonal approximation, 
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where 7( czr ) = (If) 1 / 3 , which is the same ^ Clr "> in Eq. 
(61) of Ref. 0. 

Comparing Eq. [5] with Eq. (62) of [H, we have an 
extra term due to radial periodic orbits, which is the 
single summation term in Eq. 42 of (H . 



3. variance of MK 



The variance is calculated from £(e, £7) 
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Compared with the Eq. 70 of there is the extra 
contribution from radial periodic orbits. 
The new variance equation is 
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^. numerical simulation of MK 

The ensemble averaging is averaging over an ensemble 
of /3 from a normal distribution centered at /?o = 3 x 10 6 
and with a width /3o/20. For RB, the normal distribution 
is centered at 1 with a width 0.2. 




FIG. 1: We use the same values of /3 = 3 x 10 6 , and e = 2 x 10 5 
same as in Fig. 2 of Ref. [J. The red line is the numerical 
result. The black line the old theoretical result from Eq. [6] 
The green line is the new theoretical result from Eq. [7] 



C. variance of RB 

In the old theory, the variance is Q 
£(e;£;a) = 
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In the new theory, the variance is 
Z(e;E;a) 
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Unlike the modified Kepler problem, the average of an 
ensemble of a does not converge to the average value of 
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a. Hence the correct variance should be given by 3Q 3 



E(e; E; a)p(a)da, 



which does not converge to E(e; E\ coverage)- 
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FIG. 2: The red line is the numerical result. The blue line is 
the theoretical result without ensemble averaging over a. The 
green line is theoretical fit using Eq. [9] with ensemble averag- 
ing. The black line is theoretical fit using Eq. [8]with ensemble 
averaging. We see that the ensemble averaging fits well in all 
the range, while it fits well only for E < 3000 without ensem- 
ble averaging. The variance with coherent interference is a 
little better than without. 
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FIG. 4: Rigidity of RB as a function of the energy interval 
width, e = 10 5 . 



B. rigidity of MK 



The saturation level rigidity is calculated by 
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Compared with the Eq. 69 of [l|, the extra contribu- 
tion to from radial orbits 
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III. LEVEL RIGIDITY 



A. rigidity as function of energy interval width 
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FIG. 3: Rigidity of MK as a function of the energy interval 
width, e = 2 x 10 5 . 



1. long term trend of rigidity 
For M r = 1 and (ff ) 1/3 > 1, 
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But for @ = 3 x 10 6 , this requires e » 3 x 10 6 
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2. discrepancy of rigidity 

In Fig. O there is discrepancy between numerical and 
theoretical results of rigidity. We explain it in this sec- 
tion. 

The new rigidity is that 
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C. rigidity of RB 

In the old theory, the rigidity is 
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In the new theory, the rigidity is 
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FIG. 6: The rigidity of RB. The red line is the numerical 
result. The black line the old theoretical result from Eq. 1161 
The green line the new theoretical result from Eq. |TT| The 
inset shows the rigidity change against energy width and e = 
10 5 . 



Appendix A: Derivation of correlation function, 
variance and rigidity 
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FIG. 5: We use the same values of /3 as in Ref. Q. The red 
line is the numerical result. The black line is the theoretical 
result from Eq. 1121 The green line is the theoretical result 
from Eq. [15] 



IV. CONCLUSIONS 

We presented a complete quantum mechanical deriva- 
tion of the level correlation function of a modified Kepler 
problem. With the inclusion of a class of periodic orbits 
omitted in [l[ (as well as number of other corrections), it 
is in full agreement with the periodic orbit theory. 

The account of long-range level correlations leads to 
striking properties of the level number variance and sat- 
uration level rigidity, as given by eqs. (O and (fTo) . The 
latter are the main results of this work and are in excel- 
lent agreement with the numerical calculation, as shown 
in Figs. [T] and[5J 



1. fluctuation of the level density 

The level density in the variable k has two parts 
p(k) = p (2) (k) +p (1) (k), where k = (l,p), p (2) (k) is 
a delta function at lattice points in the first quadrant 
("1Q") and half delta function at lattice points on the 
half-axes, p^-^k) is half delta function at lattice points 
on the half-axes. The level density in the variable e is 

p(e)« / d 2 kS(e - #k))p(k), (Al) 

JlQ 

where </>(k) = (2puj + 1 2 ) 3 / 2 To use the Poisson sum- 
mation formula, we make p(k) fully periodic by sym- 
metrizing the integral to the whole plane ( "WP" ) . Then 

p{e)*\( d 2 k<5( e -<//(k))p(k), (A2) 

where </>'(k) = (2|p|w + Z 2 ) 3 / 2 /3w. 

We consider p( 2 '(e) first. After applying Poisson 
summation formula, p^ 2 -* (e) becomes a summation over 

(2) 

Fourier component /w with w = 2Tr(Mg, M r ) with inte- 
ger Me and M r : 

p^\e)= ]T 4 2 >( £ ). (A3) 

M e ,M r 

A Fourier component of p*- 2 ' (e) looks like 

^(^ = I / rffc H lV 1./ "^ ( A4 ) 

where L e is the constant energy surface. For the case of 
Me = M r = 0, 

^v-f^wa- 1 (A5) 
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which is the average level density. The w = component 
gives the average level density and other components give 
the fluctuation 5p(e). The fluctuation of level density is 
given by 

Sp (2 He)=( E + E + E W(e). 

\Mg,M r ^0 M a =0,M rl tO M e ^O,M r =oJ 

(A6) 



a. case of Mg,M r > 
For every w = 2ir{Mg, M r ) with Mg, M r > 0, there are 

(2) 

four symmetric ws: 2n{±M e ,±M r ). We redefine /^(e) 
as the sum of these four terms: 
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where Mg, M r > 0. Also notice that the integral is sym- 
metric between positive and negative I or p. Using this 
symmetry, we return to the first quadrant. 
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where 
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C(ir) and S^x) are the Fresncl integrals defined as g(x) 
Jq dt ^(fi 2 ). The Fresnel integrals are almost step 
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function. We treat C+(e) = S+(e) = \, if > 1, 



which is true for large /3. 

The terms C_(e) and S'-(e) are important because 
they are almost step functions and they are the underly- 
ing origin of quantum jumps. If 



\/ T |^(M,(3we) 1 / 3 - Mgu) < 1, (A9) 

both C_(e) and SL(e) arc both C_(e) + C+(e) and 
5_(e) + S-|-(e) are 0, which makes the oscillation terms 
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not important. 

To simplify the problem, we approximate C_ (e) 
C+ (e) and 5- (e) + S+ (e) using the step function. 
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a. Quantum jumps The quantum jumps happen 
when 
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For p = 3 x 10 6 , when M e = 1 and M r = 2, e = 2.5 x 10 5 . 
When Me = 1 and M r = 1, another big jump will happen 
at e = 2 x 10 6 , which is not shown in Fig. 1 of pj. When 
M e = 1 and M r = 3, e = 0.74 x 10 5 . This jump is also 
in Fig. [5] 

b. Quantum and classical correspondence M r and 
Me correspond to winding numbers of the classical orbits 
in the angular direction and radial direction respectively. 
The orbit periods are always multiples of T m i n = —. 
When Me = 1 and M r = 1 and e = 2 x 10 6 , the system 
has reached its minimal period. 



b. case of M e = 0, M r > 
There are two symmetric terms in Eq. IA7I 
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Classically, I^Jtn m )( e ) corresponds to radial periodic or- 
bits. The rigidity before jumps comes from radial peri- 
odic orbits. 



c. case of M e > and M r = 
There are two symmetric terms in Eq 
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Classically, I^kCMo o) ( e ) corresponds to isolated angular 
periodic orbits After collecting Eq. IA121 IA141 IA151 we 
get Eq. El 
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stable isolated periodic orbits Comparing 
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f e) is far smaller than other terms. The satura- 
tion rigidity due to these stable isolated periodic orbits 
is 
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which is smaller than 6 x 10 -4 at e = 5 x 10 5 . The reason 
is that the periodic orbit of M r ^ actually denotes a 
family of continuous periodic orbits with the same M$ 
and M r . The superposition of a family of periodic orbits 
contribute the fluctuation of levels far more significantly 
than stable isolated periodic orbits. Q 
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The sin term still has no interference. Its contribution to 
rigidity is the same as the previous theory. Comparing 
with Sp^ (e), the only difference is that we have an extra 
factor 1 + \fi^- In the same way, in variance and rigidity, 

we need to add an extra factor like (l - 
corresponding terms. 

3. calculation of RB 

A similar calculation gives 
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or 4 2) = 2J (2TTM 2 k y ). lf ] = J o (0) = 1, which is the 
average spectral staircase. The two parts of is 
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2. coherent interference 



for M = Mi or M 2 . 



Now we calculate the Balian-Bloch term p^^k). Be- 
cause the semiclassical result is calculated from the as- 
sumption of fully periodic p(k), the contribution from 
5p^(e) is non-semiclassical 

On the Z-axis, 
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where w v = 2irv. 
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We can ignore this contribution to the final variance and 
rigidity. On the p-axis, 
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Sppla^) = E (3^)1/3 C ° S ^ 

A 3 (e) = i. 

Although 5pf} axia {e) is coherent with 4 2 J (Mgfi) (e), 
the coherent contribution to rigidity is small because 

(2) 
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a. coherent interference 



Bessel function J has the asymptotic form 
M x ) ~ \ —cos(x- j). 

V 7TX 4 

The only coherent interference term comes from 
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M(^) 1 / 4 ) cos 2ttM^ 



for M = Mi or M 2 . For high energy, the coherent ef- 
fect between p^> and p^ is suppressed by a factor -^74. 
The above factor is smaller than 1. This means that the 
coherent interference decreases the rigidity and variance. 



Appendix B: Ensemble averaging 

1. parametric averaging 

The ensemble average of a physical quantity can be 
mathematically defined as an integral over the distribu- 
tion function of /? or uj. The ensemble average of x is 



x(ui)p(uj)duj 



(Bl) 
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We calculate the ensemble average by sampling u> from a 
normal distribution. 



2. average spectral staircase 

After ensemble averaging, 

(Sp(e)) = 0. 

This requires that 
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From another perspective, • in cos(-) can be linearized 
to be a linear function of uj. The integration of cos(-) 



over uj gives zero. The ensemble averaging in numerical 
calculation corresponds to the diagonal approximation in 
theoretical calculation. 



3. diagonal approximation 



Here we give a new and general argument of diagonal 
approximation. Though our argument is specific to MK, 
it seems to be applicable to other problems too. There 
are many interference terms in correlation function of 



levels, such as /, 
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Mg 2 or M ri ^ M r2 . The interference between Mg 
M 02 or M ri ^ M r2 is 
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Compare with Eq. IB31 the two terms of cos of the 
above equation have the same form. The first term is 
always zero after ensemble averaging. For the second 
term, if M ri ^ M T2 and Mg 1 ^ Mg 2 , it is also zero. 
This implies that the interference between two terms with 
M ri M T2 and Mg 1 ^ Mg 2 is zero after ensemble aver- 
aging. If M ri 7^ M r . 2 and Mg 1 = Mg 2 , the form of the 
second term is also similar to Eq. IB3I So it must be zero 
after ensemble averaging. If M ri = M T2 and Mg 1 ^ Mg 2 , 
for ei ^ e 2 j the second term is zero after ensemble aver- 
aging. But it is indeed somewhat different because here 
the factor 

eT ~ (B5) 

is small. This may makes Eq. IB3I not applicable any 
more. But the main part of the second term is that 

{ttM^M T2 -nM 2 2 M ri )u; 2 _ [TiM 2 M r2 --KM 2 M ri )uj 



M ri M r2 UJ 



M ri M r2 



(B6) 

quickly oscillates with uj and becomes zero after ensemble 
averaging. 

In summary, the interference terms between different 
periodic orbits with M ri ^ M T2 or Mg t ^ Mg 2 must be 



zero. The ensemble averaging removes the interference 
between different orbits and makes the correlation func- 
tion, variance and rigidity become a sum over classical 
orbits. In this way, we give a "proof" diagonal approxi- 
mation of MK from the perspective of ensemble averag- 
ing. 



4. diagonal approximation for RB 

A typical non-diagonal term in the level correction 
function looks like 
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After replacing x — a 1 / 4 , the derivative of x in cos(-) is 
- Mi 4/^+ -^M 2 J^ > 0. The integration over a 
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gives 0. It works in the same way for the sin term with 
Mi £ M 2 . 

When Mi = M 2 , the ensemble averaging over a does 
not give 0. It is not similar to a term in 5p^ as its 
derivative of x vanishes. This term gives for rigidity. 
But its contribution to correlation function and variance 
is not zero. This needs further investigation. 

Appendix C: Scaling of variance of MK 
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FIG. 7: Variance scaled. /3 is j3 = 1, 3, 5, ■ • ■ , 19 x 10 6 . The 
black line is the theoretical result without considering p (1 K 
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